Existence of a novel metallic ferromagnetic phase in models for undoped manganites 
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The existence of a novel metal-insulator transition in the ferromagnetic state of models for un- 
doped manganites is here discussed using numerical techniques applied to the e g -orbital degenerate 
Hubbard model tightly coupled with Jahn- Teller distortions. The ground-state phase diagram is 
presented in the plane defined by the electron-phonon coupling A and Coulomb interaction u. In 
contrast to the standard one-band Hubbard model for cuprates, the metallic phase is found to exist 
for finite values of both A and u in the present e g -orbital Hubbard model even at half-filling, due to 
the Fermi-surface topology which is incompatible with the staggered orbital ordering concomitant to 
the insulating phase. Based on the present results, a possible scenario for Colossal Magneto- Resistive 
effect is discussed in undoped manganites. 

PACS numbers: 75.30.Vn, 71.30.+h, 71.10.Fd, 71.38.-k 
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One of the most interesting issues in recent research 
activities in condensed matter physics has been the eluci- 
dation of the mechanism of ColossaLMagneto- Resistance 
(CMR) phenomena in manganites BBlI Recently, it has 
been widely recognized that the CMR effect appears 
when the ground-state of manganites can be converted 
from insulating to ferromagnetic (FM) metallic by a small 
magnetic field. The stability of the FM metallic phase 
in manganites is understood due to the double-exchange 
(DE) mechanism, based on the strong Hund's rule cou- 
pling between mobile e g electrons and localized t2 g spins. 
On the other hand, the insulating phase in manganites 
occurs due to the coupling between degenerate e g elec- 
trons and Jahn- Teller (JT) distortions of MnOg octahe- 
dra, leading to various types of charge and/or orbital 
ordering, as observed experimentally.tl 

In undoped RMn03 (R=rare earth lanthanide ions), 
the parent compound of CMR manganites for most R- 
ions, the A- type antiferromagnetic (AFM) insulating 
phase appears with the,C-type ordering of the (3x 2 —r 2 )- 
and (3y 2 — r 2 )-orbitals.Q Quite recently, for R=Ho, a novel 
AFM phase called the "E-type" sraiti structure phase has 
been reported as the ground-state^Ha but this state is also 
insulating wittuthe same orbital ordering as that of the 
A-AFM phase.Ll In general, the metallic phase in man- 
ganites has been widely considered to appear only after 
holes are doped into such insulating undoped materials. 
Namely, by substituting La by alkaline earth ions such 
as Sr and Ca, holes are doped into the e g -electron band 
and, due to the DE mechanism, the FM metallic phase 
is generated. Most of the discussion in manganites has 
been concentrated on the phases induced by hole doping 
into the A- or E-AFM state with different values in the 
bandwidth of carriers. In this framework, it is implicitly 
assumed that the undoped system is always insulating. 

However, in the theoretical phase diagram for undoped 
manganites, which has been obtained by the numerical 
analysis of the DE model coupled with the JT distor- 
tions, a FM metallic phase has.-be£n observed for weak 
electron-lattice coupling region nUm This phase did not 
receive much attention since the focus of previous pa- 



pers was on the insulating phases. In fact, while it is 
an interesting possibility to have a metallic phase for 
RM11O3, its existence has not yet been confirmed even 
from the theoretical viewpoint, since the previous calcu- 
lations were carried out in small-size clusters and in addi- 
tion the Coulomb interaction was not included explicitly. 
If the FM metallic phase is confirmed to exist adjacent to 
the A- or E-AFM insulating state in the undoped limit, 
the CMR effect could in principle occur even without hole 
doping. Here, recall that at the heart of the CMR phe- 
nomena is the two-phase competition between FM metal- 
lic and insulating phases flu The insulating phase in un- 
doped manganites has been widely recognized both from 
experimental and theoretical investigations, while the ex- 
istence of the FM metallic phase has not been explored, 
in spite of its potential importance. 

In this paper, it is attempted to show the existence of 
the FM metallic phase in undoped manganites. For this 
purpose, the spinless e g -orbital Hubbard model tightly 
coupled to JT distortions is analyzed by using numerical 
techniques such as relaxation method for distortions and 
Lanczos algorithm for exact diagonalization. For small 
electron-phonon coupling or Coulomb interaction, an or- 
bital disordered phase specified with metallic character- 
istics is observed in two and three dimensions, while in 
one dimension, the orbital-ordered insulating phase is ob- 
tained except for the non-interacting case. In the phase 
diagram depicted in the plane spanned by the electron- 
phonon coupling and Coulomb interaction, the metallic 
phase is found to exist in a wide range of relevant pa- 
rameters. Then, it is concluded here that the metallic 
phase should exist even in undoped manganites, and its 
experimental confirmation by suitable chemical substi- 
tutions for RMn03 should represent a new challenge to 
experimentalists of Mn oxides. 

Models for manganites include five important ingredi- 
ents such as the kinetic term for the e g electrons, the 
Hund coupling between mobile e g electrons and local- 
ized ti % spins, the electron-lattice coupling between e g 
electrons and distortions of the Mn06 octahedra, the 
Coulomb interaction among e» electrons, and the AFM 
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coupling between neighboring t2 g spinsH Fortunately, for 
the purposes of this paper it is not necessary to consider 
all these interactions. Namely, since the FM phase is our 
focus here, e g -electron spins can be assumed to be per- 
fectly polarized from the outset due to the strong Hund 
coupling between e g electrons and <2g spins. Then, in the 
present work a spinless model including only charge and 
orbital degrees of freedom can be used. 

The Hamiltonian H studied here is, thus, given by the 
combination of three terms: i?kim -^ci-ph: and -ff e i-ei- 
The first term indicates the hopping motion of e g elec- 
trons, written as 

-ffkin = — ^ i^ T /d; r di+ar', (1) 
iarr' 

where d- la (c?ib) annihilates an e g -electron in the d x 2_ y 2 
{d 3z 2_ r 2) orbital at site i, a is the vector connecting 
nearest-neighbor sites, and t* T , is the hopping ampli- 
tude between r- and r'-orbitals along the a-direction, 
expressed as i aa = — >/3i* b = — V3t* a =3t* h =3t / 4 for a=x, 
^a=V3< b = >/3t£=3t£ b =3i/4 for a=y, and tg b =t with 
C=*ab=*ba=° for a=z > respectively. 

The second term iJ e i_ p h indicates the coupling of e g 
electrons with distortions of MnC-6 octahedra, given by 

#cl-ph = ^y^(Qli^i + Q2iT xi + Q 3 jT z i) 
i 

+ (m^iPQli + Ql + Ql), (2) 

i 

where A is the dimensionless electron-phonon coupling 
constant, Qu is the breathing-mode distortion, and Q21 
and Qsi are, respectively, the (x 2 —y 2 )- and (3z 2 — r 2 )- 
type JT-mode distortions, n\= d\ a di a + d\ h dih, r x i— 
<4 a dib+ dl h di a , and r z i= d\ a di a ~ d\ h dib- The second 
term in f/d_ p h is the usual quadratic potential for dis- 
tortions and [3 is the ratio between the spring constants 
for breathing and JT distortion mode distortions- In this 
paper, the distortions are treated adiabaticallyEJ 

The third term -ff e i-ei includes the Coulomb interac- 
tions, and it is given by 

H i_ c i = u ^2 n ia n ih , (3) 

i 

where ni T =d\ T di T and u is the on-site Coulomb interac- 
tion. Note that u=U'—J, where U 1 in the inter-orbital 
Coulomb interaction and J is the Hund's rule coupling in 
the standard notation for the multi-orbital model. In or- 
der to investigate the possibility of metal-insulator tran- 
sition in if, the model can be analyzed from many per- 
spectives, due to the combination of these three terms. 

Let us start our discussion considering the JT model, 
defined as H jT=-ffkin+-ffei-ph, obtained by dropping 
fi i_ c i. Since in this paper undoped manganites with 
one e g electron per site are analyzed, it is clear that the 
oxygen ions that are shared by adjacent MnC-6 octahe- 
dra are all active, and the concomitant distortions are 



not independent from site to site. As a consequence, the 
cooperative effect is crucially important in this case. To 
consider this issue carefully, the simplest way is to op- 
timize directly the displacement of shared oxygen ions 
by using relaxation techniques J-3 In practice, considering 
sites i and i + a, the oxygen in between is only allowed 
to move along the a-axis (i.e. buckling and rotations are 
simply neglected). Note that in one and two dimensions, 
some of the oxygen ions are not shared by neighboring 
octahedra. In this paper, for simplicity, the positions 
of such unshared oxygen ions are assumed to be fixed, 
since in the real materials those ions are not allowed to 
move freely due to the effect of compensating cations. 
In other words, the one-dimensional (ID) chain or the 
two-dimensional (2D) plane are not isolated, but they 
are embedded in a three-dimensional (3D) environment. 

In Fig. 1, the main results for Hj^ are summarized. 
For the ID chain, it is emphasized that the JT distor- 
tion occurs for infinitesimal values of A, consistent with 
Peierls instability concepts. This issue can be clearly un- 
derstood by the opening of the Peierls gap in the band 
structure, as shown in Fig. 1(a). In this case, the or- 
bital ordering with the alternation of 3x 2 — r 2 and y 2 —z 2 
orbitals, as shown in Fig. 1(b), is stabilized due to the 
cooperative JT distortions Jj On the other hand, for the 
2D case the metal-insulator transition can be observed 
at a finite value of A. To observe this, it is conve- 
nient to monitor the average value of the JT distortion 
(Q)=(l/N) J2i \JQ\\ + Qj v where N is the total number 
of sites. As shown in Fig. 1(c) using a 4x4 cluster, (Q) 
changes from zero to a finite value at A^1.2, a clear indi- 
cation of a metal-insulator transition. In this situation, 
the orbital-ordering pattern becomes bipartite as shown 
in Fig. 1(d), indicating that Q^x—I for the "A" sublattice, 
while Q2i~—q for the "B" sublattice, with Q3i=0 for all 
sites. In other words, the distortion can be expressed by 
using a single parameter q, since it occurs in a cooper- 
ative manner and all distortions are correlated with one 
another. Thus, our task can be drastically reduced in the 
numerical optimization of distortions, since it is enough 
to minimize the total energy as a function of q. Such 
a calculation has been also performed and the result is 
expressed as the dashed-curve in Fig. 1(c) for 4x4 clus- 
ter, in very good agreement with the results for the full 
optimization within the computational error-bars. 

To obtain the 4x4 cluster result of Fig. 1(c), anti- 
periodic boundary conditions are imposed to satisfy the 
closed-shell condition for the electron configuration in the 
non-interaction case. The small lattice size and use of 
special boundary conditions may raise concerns about 
the existence of the metal-insulator transition in the bulk 
limit. To clarify lattice size effects, calculations in mo- 
mentum space have been performed by imposing the 
twisted boundary condition with the Bloch phase on the 
unit cell of 2x2 or 4x4 lattice. The difference in the size 
of the unit cell appears in the sharpness of the transition, 
indicating that the period of the distortion becomes large 
at and near the transition region. As shown by solid cir- 
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FIG. 1: (a) Energy band structure of the ID chain at A=l. 
(b) Orbital ordering for the ID case, with the alternation 
of 3a; 2 — r 2 and y 2 —z 2 . (c) Average distortion as a function 
of A in the 2D plane. Solid circles and squares denote full- 
optimized results for the 4x4 cluster and in momentum space, 
respectively, while solid and dashed curves denote the results 
by using a single parameter q for the 4x4 cluster and in mo- 
mentum space, respectively. Dotted curve indicates the full- 
optimized result in momentum space using the 2x2 lattice 
unit, (d) Orbital ordering in the 2D plane, with the staggered 
pattern of Q2-type distortions, (e) Fermi-surface lines of the 
2D lattice. Solid and broken lines indicate the Fermi-surface 
line for upper and lower band, respectively, (f) Average dis- 
tortion as a function of A for the 3D lattice. Solid circles with 
curve denotes the full-optimized results in momentum space 
by using 4x4x4 unit cell, while dotted curve indicates those 
for the 2 x 2 x unit cell, (g) Orbital ordering in the 3D lattice. 



cles in Fig. 1(c), again we observe clear indications of 
a metal-insulator transition, although the critical value 
is reduced. This fact suggests that the present metal- 
insulator transition is not spurious of a finite-size cluster. 
As mentioned above, in the 2D case, the distortion has 
been expressed as a single parameter q. This is also true 
in the results for the momentum space calculation, as in- 
dicated by the solid curves in Fig. 1(c), again in good 
agreement with the full-optimized results. 

Let us consider the reason why the metallic phase can 
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u (units of t) u (units of t) 

FIG. 2: Orbital correlation function T z (q) as s function of u 
(a) in the 2D plane with 10 sites and (b) in the 3D cube with 
8 sites. 



exist even at half-filling. To clarify this point, it is quite 
useful to depict the Fermi-surface. As shown in Fig. 1(e), 
the nesting vector is (tt, 0) or (0,7r), not (tt,tt). These 
nesting vectors are not compatible with the staggered 
orbital ordering pattern (Fig. 1(d)) that is stabilized in- 
creasing A. This is one of the remarkable features of 
the multiorbital e g -electron system, which is not spe- 
cific to two dimensionality. In fact, in the results for 
the 3D shown in Fig. 1(f), we also observe the 

metal-insulator transition at a finite value of A. The dot- 
ted curves are obtained by using the 2x2x2 unit cell 
where the twisted boundary condition with the Bloch 
phase is imposed, while the solid circles are results for a 
4x4x4 unit cell. Except for the sharpness in the tran- 
sition, in both cases, a signal for a metal-insulator tran- 
sition is clearly obtained. In this case, the orbital or- 
dering pattern becomes very complicated, as shown in 
Fig. 1(g). Note that this pattern repeats periodically on 
lattice larger than 2x2x2. In the 3D case, an intrinsic in- 
compatibility between the Fermi surface and the orbital 
ordering pattern is also found. Even without invoking 
the numerical results discussed before, the qualitative ar- 
guments presented here related with the lack of nesting 
effects in Hjt at small A strongly suggests the presence 
of a metallic phase in two and three dimensions. 

Now let us turn our attentions to the effect of 
Coulomb interactions by considering the Hubbard-like 
model, i? c =-ffkm+-ffci-ci, namely dropping ff c i- P h in 
this case. This model is analyzed by exact diagonal- 
ization to measure the orbital correlation defined as 
T,(q)=(l/JV) Eije lq - (i ~ j) (r zi r,j). In the ID case (not 
shown), T z (tt) becomes dominant for infinitesimal value 
of u, indicating that the metallic phase does not exist in 
the ID systems except for the non-interacting case, as 
is well known in the standard Hubbard model at half- 
filling. However, in higher dimensions the situation is 
quite different. Figure 2(a) denotes the results for a 2D 
lattice with 10 sites, showing that the (n, tt) component 
becomes dominant for larger value of u. However, as dis- 
cussed above, the topology of the Fermi-surface is incom- 
patible with the (tt, tt) correlation, and as a consequence 
T z {ir : ir) is no longer dominant in the small-it region.Eil 
In the 3D case with 2x2x2 cube, a sudden increase of 
the (tt, 7r, tt) component is observed at u~16, since a level 
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FIG. 3: (a) E(0)—E(q) as a function of q for u~3, where E(q) 
is the ground-state energy with lattice distortions character- 
ized by a single parameter q. (b) Ground-state phase diagram 
in the (A, «)-plane. The boundary curve touches the A=0 axis 
at u«7.5. 

crossing occurs in this case. The orbital ordering pattern 
shown in Fig. 1(g) does not show (ir, tt, tt) correlation due 
to the cooperative effect. Without the cooperative dis- 
tortion, it is natural that (tt, 7T, tt) correlation becomes 
dominant, as deduced from the large- it limit. Unfortu- 
nately, since it is difficult to enlarge the lattice size for 
the exact diagonalization, it cannot be proven exactly 
that the metal-insulator transition for He occurs in the 
bulk limit. However, the results are highly suggestive 
that the staggered component in the orbital correlation 
becomes dominant in the e g -orbital Hubbard model only 
at large value of u. At small couplings, a metal should 
be stabilized. 

Although it has been quite instructive to analyze the 
effects of the electron-phonon coupling and Coulomb in- 
teraction separately, both of them are simultaneously ac- 
tive in real materials. Thus, our final task here will be 
the study of the full model £f=i?kin+-ffei-ph+-Hei-el by 
a combination of exact diagonalization for Coulomb in- 
teractions and optimization techniques for lattice distor- 
tions. It is a very hard task to perform the exact diago- 
nalization at each step during the relaxation of the lattice 
distortion, but in order to depict the ground-state phase 
diagram, it is sufficient to analyze how the total energy 
changes when lattice deformations occur. As mentioned 
in the above paragraph, fortunately in the 2D lattice the 
cooperative distortion can be expressed using a single 
parameter q. In Fig. 3(a), the total energy difference is 



presented as a function of q for several values of A at 
it=3. These results clearly indicate that the lattice dis- 
tortion does not occur for A<0.3, while the optimized 
state accompanies the lattice distortion for A>0.3, since 
the energy minimum is located at a finite value of q. 
Thus, in order to detect the metal-insulator transition 
under reasonable CPU times, it is enough to monitor the 
change of the total energy for a very small value of q. If 
it is negative, it strongly suggests the occurence of lattice 
distortions. In actual calculations, the energy difference 
is evaluated for g=0.01. Using this procedure, the phase 
diagram for the 2D 10-site cluster is shown in Fig. 3(b). 
The results clearly indicate that the metallic phase ex- 
ists, even in this case where both electron-phonon and 
Coulomb interactions are included. Note that the transi- 
tion at A=0 appears as a cross-over behavior as indicated 
above, but in the present calculations for non-zero A, the 
boundary curve seems to converge to the cross-over point 
at u«7.5, consistent with Fig. 2(a). 

Finally, the novel possibility of observing CMR effects 
in undoped manganites is here discussed. As mentioned 
in the introductory paragraph, in the theoretical phase 
diagram for the undoped limit, the FM metallic state has 
been observed adjacent to the A- or E-type AFM insu- 
lating phaseE In the present work, it has been confirmed 
that such a metallic phase is not spuriously caused by the 
small-cluster calculations, and moreover the phase does 
not disappear by adding Coulomb interactions. Thus, 
within the theoretical framework for CMR that it is based 
on the competition between FM metallic and insulating 
phases of any nature, CMR-like phenomena could occur 
even in undoped manganites, a very unexpected and chal- 
lenging theoretical prediction. 

In summary, the metal-insulator transition for un- 
doped manganites has been discussed based on the e g 
orbital Hubbard model tightly coupled to JT distortions. 
The existence of a novel FM metallic phase at zero hole 
doping has been shown, leading to the intriguing possi- 
bility of CMR effects even in RMn0 3 . 
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